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Abstract. We present the complete set of ISOPHOT observations of 3CR radio galaxies and quasars, which are contained 
in the ISO Data Archive, providing 75 mid- and far-infrared spectral energy distributions (SEDs) between 5 and 200 nm. For 
28 sources they are supplemented with MAMBO 1.2 mm observations and for 15 sources with new submillimetre data from 
the SCUBA archive. The sample includes flat and steep spectrum quasars, broad and narrow line radio galaxies, as well as 
Fanaroff-Riley FR 1 and FR 2 types. The SED shapes exhibit a diversity in the infrared (IR), ranging from a smooth dominating 
synchrotron component in flat spectrum sources to a thermal dust bump around 60-100 /im in steep spectrum sources. The 
detection rate of a thermal bump in more than 50% of the cases suggests that dust emission may be a general phenomenon 
in these sources. We check the orientation-dependent unified scheme, in which the powerful FR2 narrow line galaxies are 
quasars viewed at high inclination, so that their nuclei are hidden behind a dust torus intercepting the optical-ultraviolet AGN 
■»^ , radiation and reemitting it in the infrared. We find that (1) both the quasars and the galaxies show a high mid- to far-infrared 

rG ' luminosity ratio typical for powerful AGNs and (2) - when matched in 178 MHz luminosity - both show the same ratio of 

M-(| isotropic far-infrared to isotropic 178 MHz lobe power. Therefore, from our large sample investigated here we find strong 

Q . evidence for the orientation-dependent unification of the powerful FR 2 galaxies with the quasars. The distribution of the dust- 

i—{ ' to-lobe luminosity ratio shows a dispersion which we suggest to be most likely due to the additional influence of evolution 

jy^ , and environment superposed on the orientation-dependent unification. We discuss our data also in the frame of the receding 

C^ ■ torus model. At the high 178 MHz luminosities of our sources above 10'^'^"'' W/Hz we do not find any support for this model 

in its original formulation and therefore we propose a refinement: The scale height of the torus might not be independent of 
luminosity, rather it may increase at high luminosities due to the impact of supernovae from starbursts accompanying the AGN 
J^ ' phenomena. 

^ Key words. Galaxies: fundamental parameters - Galaxies: photometry - Quasars: general - Infrared: galaxies 

1 . Introduction (3) the aspect angle. To disentangle these effects is one of the 

major challenges of the current AGN research, as reviewed by 
1. 1. Unified schemes eg urry & Padovani (1995). 

While observationally at least a dozen classes of active galactic The Revised Third Cambridge Catalogue of radio galax- 

nuclei (AGN) are discernible, on the conceptional side consen- ies and quasars at 178 MHz by Laing et al. (1983, abbrevi- 

sus is growing that in most, if not all, of these different objects ated 3CRR) and - applying softer criteria - by Spinrad et al. 

the accretion onto a massive black hole (Mbh > 10^ M©) pro- (1985, for short 3CR) contains 178 and 298 sources, respec- 

vides the main source of energy. It seems natural to explain tively, which fall into several classes of radio-loud AGN: 

their apparent differences by observational effects such as (1) 

the circumealactic environment, (2) the evolutionary phase, or , , „ ... , , • , 

1) flat spectrum quasars, which are almost core dominated, 



o 
o 



Senrfr#pn>ifre9««r.. to: Martin Haas (haas@astro.rub.de) 2) steep spectrum broad Hne radio galaxies (BLRGs) and 

* Based on observations with the Infrared Space Observatory ISO, quasars (here collectively called quasars), with either com- 

an ESA project with instruments funded by ESA Member States (es- pact or lobe dominated morphology, 

pecially the PI countries: France, Germany, the Netherlands and the 3) powerful narrow line radio galaxies (NLRGs) of Fanaroff- 

United Kingdom) and with the participation of IS AS and NASA. Riley type FR2 with - often asymmetric - jets ending in 
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edge-brightened double lobes extending far beyond the host 
galaxy, 
4) low power narrow line radio galaxies of Fanaroff-Riley 
type FR 1, with double jets feeding edge-darkened radio 
lobes often contained within the host galaxy. 

The diversity and the brightness of the sources makes the 3CR 
catalogue a well suited database for studying the relationship 
between these AGN classes. The aspect angle relation between 
flat and steep spectrum quasars, with the flat spectrum sources 
being steep spectrum sources seen pole-on (as proposed by Orr 
& Browne 1982), is now widely accepted. The relation between 
steep spectrum quasars and FR2 radio galaxies and that be- 
tween FR 1 and FR 2 sources is still a matter of debate. 

While the quasars show the optical signatures of a powerful 
AGN (like broad emission lines and high excitation and opti- 
cal luminosity), the galaxies do not. In an early attempt, phys- 
ically different central engines were proposed as reviewed by 
Begelman, Blandford & Rees (1984): quasars have a high mass 
accretion rate but a low black hole spin, and galaxies have a low 
mass accretion rate but a high black hole spin. Both engines 
should be able to create the radio jets. Although a viable expla- 
nation, from the conceptional view point it appears somewhat 
astonishing that two intrinsically different AGN types should 
lead to the same radio jet and lobe phenomena. Alternatively, in 
the picture of orientation-dependent unification, galaxies could 
be quasars in which the AGN is hidden behind a dust torus 
seen "edge-on", as proposed by Barthel (1989). Further debate 
arises from the fact that galaxies are distributed over the en- 
tire redshift range of the 3CR catalogue from the nearby lo- 
cal universe at z « 0.01 up to z « 2, but the quasars are prefer- 
entially found at larger distances at z J5 0.3 and 178 MHz lu- 
minosities above lO^*'-^ W/Hz (e.g. Singal 1996, Wiflott et al. 
2000). The increase of the quasar fraction with distance and lu- 
minosity may find an explanation by the receding torus model 
(Lawrence 1991, Hill et al. 1996): due to dust sublimation the 
inner wall of the torus recedes with increasing luminosity of 
the central engine and - under the assumption that the torus 
scale height is independent of luminosity - the cone opening 
angle increases, leading to a higher chance to see the broad line 
region. However, it is still under debate how far this elegant 
model is valid for the powerful radio galaxies with 178 MHz 
luminosities above \(fi^-^ W/Hz. Alternatively, also physical 
or cosmic evolution could play a role: in this picture galaxies 
or a subset of them could be old quasars with a starving black 
hole, but still with luminous radio phenomena, possibly due to 
a favourable circumgalactic environment. Support for an envi- 
ronmental or evolutionary relationship comes from the fact that 
many FR2s and most FR Is fall into the low ionisation emis- 
sion line category (Hine & Longair 1979), with little evidence 
for thick obscuring dust tori on HST images (Chiaberge et al. 
1999, 2000). 

In order to test the orientation-dependent unification be- 
tween FR 2 galaxies and quasars, henceforth simply called uni- 
fication, the task is to reveal a hidden powerful AGN in a (pow- 
erful) galaxy. Clear evidence for an AGN hidden behind a dust 
torus can be obtained by spectropolarimetry, when broad lines 
show up as light scattered in the bipolar cones with the polar- 



isation angle being perpendicular to the polar axis (Antonucci 
& Miller 1985). Despite a number of cases with broad lines re- 
vealed by spectropolarimetry, in numerous cases huge extinc- 
tion or dust lanes extended over kpc could have prevented the 
view into the central region via a "suitable mirror" (Antonucci 
& Barvainis 1990, Hifl et al. 1996, Cohen et al. 1999). Note 
that the optical signatures of dust extinction only indicate that 
something could be hidden, but it remains to be shown that 
there is actually a powerful energy source behind the dust lanes. 
A robust check of the unification is to look for the mid- and 
far-infrared reemission of the absorbed light from the AGN. A 
great advantage is that at wavelengths A^25 /im the IR emis- 
sion is largely optically thin, hence isotropic and independent 
of the aspect angle. 

1.2. Infrared Observations 

Observations of the 3CR sample with IRAS did not allow a 
substantial conclusion about the unification (Heckman et al. 
1992, Heckman et al. 1994, Hes, Barthel & Hoekstra 1995, 
Hoekstra, Barthel & Hes 1997). The wavelength coverage and 
sensitivity of the IRAS observations were not sufficient to dis- 
criminate between dust and synchrotron emission in the in- 
dividual sources, hence to correct for the contribution by a 
beamed synchrotron component to the mid- and far-infrared 
spectra of the sources. 

With the Infrared Space Observatory ISO (Kessler et al. 
1996) several small 3CR sub-samples were successfully ob- 
served: On the basis of a few ISOPHOT-MAMBO detections 
Haas et al. (1998) demonstrated that the IR emission of quasars 
and galaxies is a mixture of thermal and synchrotron radi- 
ation, dominated by thermal emission in radio-quiet quasars 
and synchrotron emission in flat spectrum radio-loud ones. 
Remarkably, steep spectrum quasars as well as galaxies showed 
powerful infrared dust emission. 

In order to test the unification, a dedicated strategy com- 
pared galaxy - quasar pairs, which match in isotropic 178 MHz 
radio luminosity as well as redshift (in order to minimize pos- 
sible cosmological effects). However, the four FR2 galaxy - 
quasar pairs studied by van Bemmel, Barthel & deGrauuw 
(2000) yielded a lower detection rate of thermal dust emis- 
sion and a lower dust luminosity for the galaxies than for the 
quasars, contrary to the expectations from the unification. On 
the other hand, ISOPHOT-MAMBO observations of ten FR2 
galaxy - quasar pairs with high radio power by Meisenheimer 
et al. (2001) yielded a balanced detection statistics and a sim- 
ilar dust luminosity for the quasars and galaxies, as predicted 
by the unification. Furthermore, the extrapolation of the syn- 
chrotron spectra of the quasars from the cm and mm regime 
into the far- and mid-IR has on average a more than ten times 
higher level than that for the radio galaxies. This is also con- 
sistent with the unified schemes, which predict that the jet axis 
of quasars is oriented more toward the line-of-sight, resulting 
in a stronger beamed component. Further ISOPHOT-MAMBO 
studies of small 3CR samples confirmed that many galaxies 
as well as quasars produce luminous dust emission (Andreani 
et al. 2002), and that the synchrotron component in quasars 
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is higher than in galaxies (van Bemmel & Bertoldi 2001). A 
closer quantitative look at the dust-to-radio luminosity ratio, 
however, reveals also differences between galaxies and quasars 
which suggest the additional role of evolution and environment 
superposed on the unification (Meisenheimer et al. 2001, Haas 
2001). 

Due to the limited detection of individual 3CR sources by 
IRAS, no definite conclusions about the FR1-FR2 relation 
could be drawn; to our knowledge neither ISO results on that 
topic have been reported so far. 

1.3. This paper 

In order to overcome the small number statistics and derive re- 
sults from a broader statistical basis, we present here the re- 
maining unpublished ISOPHOT SEDs of the 3CR sources in 
the ISO archive. We discuss the statistical properties of the 
complete mid- to far-infrared data base of 75 galaxies and 
quasars from the 3CR catalogue observed by ISOPHOT with 
regard to unification, evolution and environment as well as the 
applicability of the receding torus model. We focus here on 
the unification; the discussion of the FR1-FR2 relation will 
follow in a subsequent paper. Since common selection criteria 
for the sources were low cirrus foreground and good visibil- 
ity to the satellite, this sample provides a fairly representative 
sub-sample of the 3CR catalogue, with a redshift range from 
nearby to distant (z ^ 2) objects. Throughout this paper we use 
blackbodies modified with a dust emissivity index of /3 = 2. We 
adopt a A cosmology with Hq = 71 kms^^ Mpc~^, ^matter = 
0.27 and U^ = 0.73 (Bennett et al. 2003). 



2. The data 

The ISO Data Archive (Kessler et al. 2000) contains 5 to 200 
fim photometry for 75 3CR sources obtained with ISOPHOT 
(Lemke et al. 1996), the photometer on board ISO. About half 
of them have also been observed between 5 and 15 /xm with 
ISOCAM. Our sample does not contain those 3CR sources, 
for which only ISOCAM mid-infrared data are available, with- 
out any far-infrared complement by ISOPHOT data or IRAS 
detections. The full set of 85 ISOCAM observations is pub- 
lished by Siebenmorgen et al. (2004). Our sample does not 
include the two FR 1 sources 3C71 (NGC 1068) and 3C231 
(M 82), which do not belong to the classical FR Is with ellipti- 
cal hosts, and which are nearby and therefore too extended to 
be observed in ISOPHOT photometry modes. In order to char- 
acterize the (sub)-mm part of the SEDs for the full sample, we 
complemented the available literature data: For 15 sources we 
found submm photometry in the JCMT-SCUBA archive, and 
performed additional mm observations for 28 sources with the 
Max-Planck Millimeter bolometer array MAMBO (Kreysa et 
al. 1998) at the IRAM 30 m telescope. 

2.1. ISOPHOT 

The observing modes (Laureijs et al. 2002) comprise 
chopped measurements, including those with the spectrome- 



ter ISOPHOT-S, from which we determined broad band fluxes, 
and small maps. 

We have reduced the data using the ISOPHOT Interactive 
Analysis tool (PIA' VIO), together with the calibration data set 
V7.0. Using the latest versions of the data reduction tools, there 
are no systematic photometric offsets with regard to PIA V7 
and V8 used for most of the results published earlier by various 
authors. The accuracy of the absolute photometric calibration 
depends mainly on systematic errors as described in Laureijs & 
Klaas (1999), and it is currently known to be better than 30% 
for faint sources. The relative filter-to-filter calibration agrees 
within 15% (Klaas et al. 2002). 

2.2. MAMBO 

With the MAMBO array at the IRAM 30-m telescope on Pico 
Veleta, Spain, during the pooled observation campaigns be- 
tween January and March 2003, we obtained 1.2 mm (250 
GHz) continuum flux densities for 28 of our sources. We used 
the standard on-off photometry observing mode, chopping be- 
tween the target and sky at 2 Hz, and nodding the telescope 
every 10 s. The total on sky integration time varied between 
2.5 and 23 minutes per source, and each object was observed 
at least twice, on different days. The absolute flux calibration 
was established by observations of Mars and Uranus, resulting 
in a flux calibration uncertainty of approx. 20%. The data were 
reduced using the MOPSI software package. 

2.3. SCUBA 

Yet unpublished submm continuum data obtained with SCUBA 
were retrieved from the JCMT archive for 15 sources. The ob- 
servations were performed in photometry mode, and for one 
source (3C 274 = M 87) as jiggle maps. We have reduced them 
with the SCUBA User Reduction Facility (SURF), including 
identification of noisy bolometer pixels and removal of sky 
noise. The atmospheric transmission was determined from sky- 
dips and water radiometer data at the CSO. IRCh- 10206 and 
HL-Tau served as standard calibrators. The integration time per 
source was 10-20 min. In case of very good atmospheric trans- 
mission (tcso ~ 0.05) we also included the 450 /^m data. The 
photometric accuracy at 450 /zm and 850 /im is about 30% and 
25%, respectively. 

2.4. Literature data 

ISO photometry of 49 3CR sources has so far been pub- 
lished by various authors, as labelled in col. "publ" of TableQ] 
The SED plots shown here include 4.5 to 14.3 /im data with 
ISOCAM (Cesarsky et al. 1996) for 33 sources, which are 
published by Siebenmorgen et al. (2004) and for a subset by 
Freudling et al. (2003). 

Additional literature data were retrieved from the NED and 
SIMBAD, the 2MASS archive and numerous papers: Lilly & 



' PIA is a joint development by the ESA Astrophysics Division and 
the ISOPHOT consortium led by the MPI fiir Astronomic, Heidelberg. 
Contributing institutes are DIAS, RAL, AIP, MPIK, and MPIA. 
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Fig. 1. Spectral energy distributions of those quasars, BLRGs and FR2 NLRGs which show a clear thermal bump and are used 
in the discussion of the unification in Sect. l4.2l The measurement errors are of the size of the symbols. The wavelength and 
frequency ranges are as observed and not corrected to the rest frame of the objects. The MIR (10-40/im) and FIR (40-1000/im) 
wavelength ranges, as well as the temperature of the modified blackbody (solid line curve) fitted to the thermal FIR bump are 
shown in the object's restframe. Where the data coverage is sparse, the straight lines indicate the "envelopes" used for computing 
the restframe luminosities. 
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Fig. 2. Spectral energy distributions of the FR 1 sources. The meaning of the lines and labels is as in Fig.[T] Where the ISOPHOT 
data coverage is sparse, IRAS flux values from the NED are shown, if available (e.g. for 3C 293). 



Longair (1984), Spinrad et al. (1985), Laing et al. (1983), Lilly 
et al. (1985), Neugebauer et al. (1987), Simpson et al. (2000) 
and de Vaucouleurs et al. (1991) at optical and NIR wave- 
lengths, Chini et al. (1989a,b), Steppe et al. (1988, 1992, 1993), 
Robson et al. (1998), Hughes et al. (1993), Best et al. (1998), 
van Bemmel et al. (2000), PoUetta et al. (2000), Archibald et al. 
(2001), Willott et al. (2002), Andreani et al. (2002) and Stevens 
et al. (1998) at mm wavelengths, and at radio wavelengths from 
Jodrell-Bank Web pages (at http://www.jb.man.ac.Uk/atlas/i, 
Kiihr et al. (1981), Keflermann et al. (1989) and Akujor et al. 
(1994). 



3. Results 

3.1. General properties 

The IR, submm and mm fluxes of the 3CR sources are listed in 
Tabled ISOPHOT achieved a detection rate of about 70% in 
the MIR (< 25 ^im) and 60% in the FIR (> 60 ^m). Even at the 
long wavelengths between 120 and 200 yum the detection rate 
is about 20%. The fluxes of the brightest non-variable sources, 
also detected by IRAS (e.g. 3C 272.1, 3C405, 3C459), agree 
within 15%. 

14 out of 28 sources observed with MAMBO at 1.2 mm 
were detected, and 3-0- upper flux limits are provided for the 
others. From the 15 sources in the SCUBA archive eight were 
detected at 850 /xm and three at 450 fim. 

Supplementing our data with further photometry from the 
literature as listed in section l24l we display the observed spec- 
tral energy distributions (SEDs) in Figs.n|2land|3] whereby 
each Figure contains a subset of sources suitably chosen with 
regard to the illustration of the discussion items. Fig. ^ shows 
only those quasars, BLRGs and FR 2 NLRGs which exhibit a 
clear thermal bump and are used for investigating the unifica- 
tion. Fig.|2]shows the FR 1 SEDs. Finally, Fig.|3]shows the re- 
maining SEDs which are either those of flat spectrum sources 
or which exhibit poorly constrained thermal bumps due to up- 
per limits only - nevertheless illustrating the high IR luminos- 
ity limits. Note that the SEDs in Figs.[2and|2]run smoothly 
and the photometric uncertainties represented by the size of the 
symbols are small. 



The most remarkable feature of the SEDs in Fig.[2is the 
steady increase in flux density from near- to far-IR wavelengths 
with a peak typically at 60-100 /im, followed by a steep drop at 
(sub)-mm wavelengths in many cases. Where it could be mea- 
sured, in particular for 3C 48, 3C 298, 3C 3 1 8, 3C 405, this drop 
follows exactly the Rayleigh-Jeans tail of the modified black- 
body (with dust emissivity index /3 = 2). This strongly suggests 
a thermal nature of the FIR emission. 

In order to discuss the SEDs in the framework of ther- 
mal emission, we have eyeball fitted a modified blackbody to 
the peak of the flux distribution, as shown in Figs.[2 and |2] 
The peak blackbody is mainly determined by the maximum 
in the FIR and, if available, also by the sub-millimetre points. 
However, this greybody fit does not always match the (sub)- 
millimetre measurements, since they may often be dominated 
by the synchrotron emission (for example in 3C 109 or 3C234 
in Fig.[3. Note that our "peak blackbody" corresponds in a i/F^ 
versus ly plot to the "break" of the thermal spectrum, where 
the slope equals unity. Its temperature Tp (corrected for red- 
shift) lies between ~30 and lOOK, with higher Tp values pref- 
erentially found in the most luminous sources at high redshift 
- probably a natural consequence for a flux limited sample like 
ours. It is clear from Fig.^that also dust which is colder than 
the peak blackbody may be present, but if so, then at lower flux 
levels. Also, as being evident from the MIR data points above 
the peak blackbody curve, warmer dust components up to the 
dust sublimation temperature of about 1500 K are present, but 
for lucidity of the SED plots they are not shown explicitly here. 
Note that for the cases of upper limits and poorly determined 
dust bumps shown in Fig.|3]we tentatively adopt an average Tp 
'^ 40 K determined from the good cases shown in Fig. [2 if not 
constrained sufficiently. 

The restframe thermal mid- (10-40 fim) and far-IR (40- 
1000 /im) luminosities are computed by integrating the SED 
"envelopes". These curves are either sufficiently sampled by 
observations, as e.g. in the case of 3C 48, or they are partly sam- 
pled between 10 and 100 fim and follow the Rayleigh-Jeans tail 
of the peak blackbody between 100 and 1000 /im, as e.g. in the 
case of 3C 47. For a few sources in Fig.[T](3C 33, 3C 49, 3C 245 
and 3C 284) the SEDs are sparsely sampled in the MIR-FIR 
regime; in these cases they are interpolated linearly in the log- 
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log SED plots between NIR and FIR flux points as indicated by 
the solid lines. Since the short wavelength part of the integra- 
tion interval dominates the luminosity calculation, for poorly 
sampled SEDs the MIR luminosity may be overestimated by up 
to ^50%, but the FIR luminosity is not significantly (<15%) 
affected by the choice of Tp. 

For all sources Table |2] lists the derived parameters like 
luminosities, dust temperatures and masses, and other values 
used in the discussion. 

3.2. Notes on some individual sources 

3C274 (M87): The SED of this FR 1 source clearly exhibits 
as three basic features (1) the optical-NIR bump from the host 
galaxy, (2) the thermal dust bump in the mid- and far-IR, and 
(3) the radio synchrotron component. Tabled lists the total 
fluxes at 450 and 850/im, nevertheless we were able to resolve 
the one-sided jet and the core in the SCUBA jiggle maps. The 
core fluxes assuming a point source with FWHM 8" and 15" 
are 0.391 ±0. 10 Jy and 1.085±0.10 Jy at 450 and 850 ^im, re- 
spectively. 

3C268.4: Our FIR fluxes at 90 and 170 nm are similai- 
to those reported by Andreani et al. (2002), but they are sur- 
prisingly high for a quasar at z « 1 .4. Therefore we checked 
for the presence of other nearby galaxies. In fact, the galaxy 
NGC4138 lies less than 200" apart and has an extent of more 
than 100". In the FIR it is also bright with F(170/im) > 6 Jy, 
as we found with the help of the ISOPHOT Serendipity Survey 
(Stickel et al. 2000). Hence the large ISOPHOT beams of 45" 
and 180", respectively, have seen mainly the contamination by 
NGC 4138. Therefore, we do not consider 3C 268.4 any further 

3C313: The FIR flux reported by Andreani et al. (2002) is 
due to a misidentification, actually the ISO observations refer 
to IRAS F15086H-0801 (^ 3C 313). 

3C368: The source was observed three times in the MIR 
5-15 nm range, once with ISOPHOT on Oct. 18, 1996 (at 4.8 
and 7.3 /xm) and twice at 12/xm with ISOCAM on Feb. 27, 
1996, and on Oct. 27, 1997. The detections on Oct. 18, 1996 
exceed by far the 3ct upper limits from Feb. 27, 1996 and Oct. 
27, 1997. The ISOPHOT chopped series show a perfectly os- 
cillating on-off pattern leaving no doubt about the reality of 
the source detection. Since the explanation of strong MIR vari- 
ability of 3C368 seems unlikely to us, we suggest that the 
"MIR flux excess" during the observations on Oct. 18, 1996 
is due to the passage of an asteroid (or comet), seen in the 
ISOPHOT aperture. However, no so far known asteroid has 
crossed the ISOPHOT observations (private communication by 
Thomas Miiller). Nevertheless, an asteroid (or comet) seems to 
be the best explanation for the spurious 3C 368 detection and 
non-detections. Furthermore, at a temperature of 230 K as indi- 
cated by the MIR fluxes, an asteroid superposed on the 3C 368 
field of view does not affect the FIR fluxes of 3C 368 at A > 
60 ^m. Since 3C 368 was also detected with SCUBA at 450 /.tm 
(Archibald et al. 2001), we ascribe the FIR-submm dust bump 
to 3C 368 and use the ISOCAM upper limit for MIR luminos- 
ity estimates. Nevertheless, we do not include 3C368 in the 
discussion here. 



3C 390.3: This broad line radio galaxy was also detected by 
IRAS, but both the ISOPHOT and the ISOCAM MIR data do 
not confirm the exceptionally high 12 and 25 /im IRAS fluxes 
(Miley et al. 1984), thus making the mid- to far-IR colours of 
3C 390.3 more similar to those of other BLRGs. The high IRAS 
fluxes could be contaminated by an asteroid or comet. 

3C 454.3: This is one of the few examples of a. flat spec- 
trum quasar with a strong FIR dust bump sticking out of the 
bright synchrotron continuum. Other such examples are 3C 273 
(Fig.|3}, 3C279 (Haas et al. 1998) and PG 1302-102 (Haas et 
al. 2000). We checked, whether the bump in 3C 454.3 could 
be due to another galaxy, but the only candidate on DSS and 
2MASS images is an unresolved object about 15" NNE with 
a K-star like optical-NIR SED. Therefore we conclude that the 
bright FIR dust bump is in fact caused by 3C 454.3. 

4. Discussion 

Figs.[T] and 121 demonstrate that at least half (> 40/75) of the 
3CR SEDs display a mid- to far-infrared bump due to thermal 
emission by dust. Even in the remaining cases the upper lim- 
its allow for considerable dust emission (Fig.|3] Table|2j. In the 
following we consider the different classes (e.g. flat and steep 
spectrum quasars, BLRGs, FR 2 NLRGs and FR Is), and inves- 
tigate how the dust emission compares to the properties at other 
wavelengths. First, we discuss the SED shapes for represen- 
tative class members. Then we study the statistical properties 
of the FR 2 NLRGs, quasars and BLRGs with respect to the 
unification, whereby we also consider other effects like those 
predicted by the receding torus model as well as evolution and 
environment. 

4.1. SEDsiiapes 

The SEDs shown in Figs.^andlJlexhibit a striking diversity of 
shapes. Figure@]shows five representative examples among the 
different AGN classes, i.e. radio and optical types: 

1) Flat spectrum core dominated quasars ("FSQs") like 
3C 454.3: The synchrotron spectrum extends at a nearly con- 
stant flux level from the cm to the mm range, and dechnes at 
shorter wavelengths. In the FIR, the SEDs may show a thermal 
dust bump, as in the powerful ESQ 3C 454.3, or not, as in the 
less luminous BL Lac object 3C 371 (Fig.|3}- In the context of 
unified schemes these flat spectrum sources are viewed nearly 
pole-on as proposed by Orr & Browne (1982). 

2) Steep spectrum quasars like 3C 35 1 : The synchrotron 
spectrum drops from the cm to the mm wavelengths, so that 
the thermal FIR dust emission bump becomes clearly visible. 
The optical spectra are quite blue, typical for an AGN outshin- 
ing the host galaxy. Also, the SEDs show a power-law like flux 
increase, from the NIR over the MIR and peaking in the FIR, as 
seen also in optically identified PG quasars (Haas et al. 2003). 
With respect to the optical flux, the FIR dust bump may be 
quite pronounced with a steep slope from the NIR to the FIR 
as in 3C48 at Z'--^0.3 and 3C318 at z^l.5 or more shallow as 
in 3C249.1 at Z'-0.3 and 3C298 at z-1.5 (Fig.O. We do not 
see any difference between the SEDs of FR 2, lobe dominated 
and CSS types. 
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Fig. 4. SED shapes of 3CR sources as representative examples 
of the different optical and radio classes, having good spectral 
data coverage. The symbols are as in Fig.^ 



3) Broad line radio galaxies (BLRG) like 3C33.1: The in- 
frared and radio spectra are similar to those of the quasars, 
showing in the IR also the powerlaw rise between 5 and 60 /im. 
In contrast, the optical spectrum is redder than that of quasars, 
probably due to a stronger contribution of the host galaxy rela- 
tive to that of the AGN, and possible obscuration of the intrin- 
sically blue central spectrum. 



4) Narrow line radio galaxies (NLRG) of class FR 2 like 
3C405: The far-infrared and radio spectra are similar to those 
of quasars and BLRGs. Relative to the FIR, however, the emis- 
sion in the optical and also in the NIR is weaker This sup- 
pression at short wavelengths shows some range: while for 
3C 405 the optical-NIR SEDs are only "moderately" reddened, 
the optical to IR slope is very steep for both 3C 65 and 3C 234 
(Fig.Q. In the unified scheme the NLRGs of class FR2 rep- 
resent edge-on quasars. 3C 234 shows polarised broad H^ but 
narrow H^ lines (Antonucci 1984); hence it might be seen at 
the grazing angle to the edge of the dust torus. Some NLRGs 
of class FR2 like 3C 321 and 3C459 (Fig.[lli show a strong FIR 
bump and only a moderate MIR contribution, more reminiscent 
of the SEDs of starburst galaxies. 

5) Narrow line and low-excitation radio galaxies (LERG) 
of class FR 1 like 3C 272.1, showing three basic SED features: 

1) in the optical-NIR a bright bump due to the elliptical hosts, 
represented by a blackbody of T w 3500 K peaking at about 
1 /^m, 

2) in the cm-mm wavelength range a steep synchrotron spec- 
trum, and 

3) in the MIR-FIR a well discernible dust emission bump. 

The examples shown in Fig. |3 illustrate the basic SED fea- 
tures. In the optical they arise from the host galaxy and the 
blue AGN continuum, in some cases reddened by extinction 
toward the central region. In the IR and radio ranges the emis- 
sion from dust and from the jet and lobes can well be discerned. 
Remarkably, the SEDs exhibit also a diversity of shapes even 
within each AGN class, for example in the strength of the FIR 
dust bump or the slope between the NIR and FIR fluxes. For 
the optically selected PG quasars, which are thought to be seen 
nearly pole-on, such differences in the IR emission could be in- 
terpreted as evolution of the dust distribution around the heat- 
ing sources (Haas et al. 2003). Although the diversity of SED 
shapes within each 3CR AGN class may also indicate such in- 
trinsic differences, the orientation with respect to the line of 
sight is not yet known for our 3CR sample and it is still un- 
der debate, how far the observed near- and mid-IR emission is 
influenced by different viewing angles. Therefore, we will re- 
strict our investigation to the orientation problem. Furthermore, 
average SEDs for each AGN class could be computed, but they 
will smear out the details recognisable in the individual SEDs. 
Therefore, instead of using averages, in the next subsection we 
prefer to analyse the properties from the individual SEDs. 

4.2. On the FR2 unification 

In this section we consider the following two basic classes: 
i) the steep spectrum quasars and the BLRGs, henceforth for 
short denoted as quasars, and ii) the FR 2 NLRGs, henceforth 
denoted as galaxies. The strategy to check the unification in- 
cludes two steps: (1) to show that both the quasars and the 
galaxies exhibit a high mid- to far-IR luminosity ratio typical 
for AGNs, and (2) that the isotropic FIR-to-radio luminosity 
ratio is the same for quasars and galaxies at matched isotropic 
178 MHz radio power. We consider only those sources shown 
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in Fig.Q] which have sufficiently well sampled SEDs. The sam- 
ple investigated here consists of 17 galaxies and 18 quasars, 
hence it can be considered as balanced. The set of sources 
shown in Fig. |3] and not included here comprises 20 galaxies 
and 15 - partly flat spectrum - quasars, and we will see below 
that their exclusion does not introduce a strong bias. 

4.2.1 . AGN typical dust emission 

The thermal 10-1000 fim infrared luminosities of our sample 
span four orders of magnitude, frommoderate 10^° Lq over lu- 
minous 10^^ Lq and ultraluminous lO^'^ Lq to hyperluminous 
lO^'^-lO^'* Lq objects (Tableland abscissa of Fig.|5}- The most 
luminous objects lie at the highest redshift, as is typical for flux 
limited samples like the SCR catalogue. 

The dust masses derived from the peak blackbodies shown 
in Fig.^ range between 10^ and 10^ Mq (Table|2}, which is 
comparable to those found in PG quasars (Haas et al. 2003) 
and in ULIRGs (Klaas et al. 2001). Although the derived dust 
masses are possibly uncertain by a factor of ten, they do still 
provide a rough estimate of interstellar material associated with 
these objects. 

In order to compare the 3CR sources with other samples 
of known properties, we have decomposed the IR emission 
into the MIR (10-40 A*m) and FIR (40-1 000 ^m) luminosities 
(Table|2li. For the quasars as well as for the galaxies the mid- 
to far-infrared luminosity ratio ranges between 0.4 and 5 as 
shown in Fig.|5] This is comparable to what is found for opti- 
cally selected PG quasars (Haas et al. 2003). In contrast, those 
ultraluminous infrared galaxies (ULIRGs), in which the dust 
emission is powered mainly by starbursts and not by an AGN, 
populate the Lmir/Lfir range from 0.15 to 0.5 (Klaas et al. 
2001). This is significantly lower than that of both the quasars 
and the galaxies. Therefore, the high mid- to far-infrared lumi- 
nosity ratio above 0.5 together with the high dust luminosity 
(Lift ^10^^ Lq) provides evidence for the presence of a pow- 
erful AGN in the quasars and the galaxies as well. Their AGN 
may be accompanied by starbursts, but for our purpose here the 
presence of a powerful buried AGN in the galaxies is relevant. 
It argues in favour of the unification. It should be noted here, 
that the orientation-dependent unification picture may not ap- 
ply for the FR 1 and the low -power FR 2 sources, which have 
relatively low dust luminosities and masses. 

A closer look on the mid-to-far-infrared luminosity ratio 
shows that the quasars lie above unity, while many of the 
galaxies lie below unity (Fig.|5jl. Such a difference between 
quasars and galaxies is consistent with unified models, if the 
proposed dust torus becomes optically thick at MIR wave- 
lengths (Pier & KroUk 1992, 1993, van Bemmel & DuUemond 
2003). In this case, if seen edge-on, our view towards the hot 
MIR emitting dust at the inner walls of the torus will suffer 
from extinction even at MIR wavelengths by the outer parts of 
the torus or by more extended dust. The actual degree of MIR 
opacity is not yet known. A value of Amir. ^ 1-2 is needed 
to explain the lower Lmir/Lfir ratios for the galaxies. Using 
standard extinction curves this value corresponds to Ay ^ 100- 
200, hence not exceptional with regard to that found in ULIRGs 
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Fig. 5. Mid- to far-infrared dust luminosity ratio versus 10- 
1000 /im luminosity. The rectangular box shows the range of 
local starburst ULIRGs (from Klaas et al. 2001), each plotted 
with a "H-" sign, and the dotted line indicates the upper end of 
their Lmir / Lfir distribution. The only 3CR source falling 
below this line is 3C459. 



(e.g. Haas et al. 2001). Therefore we do not yet see any con- 
flict with the unification and further investigations with more 
detailed observations as well as models may clarify that issue. 

4.2.2. Radio-Infrared-Comparison 

So far we have found evidence for a powerful AGN in both 
the quasars and the galaxies. In a strict sense, however, the 
concept of unification requires that for an object drawn from 
the parent population any isotropic emission remains the same 
while rotating the viewing angle to its axis. Thus, for an ideal 
sample of parent objects the emitted isotropic FIR dust power 
should be the same for objects of identical isotropic lobe power 
Therefore, following the ideas by Meisenheimer et al. (2001), 
we consider Rdi, the ratio of P^, at FIR wavelength 70 /im (= 
4.3 THz) and P^ at radio frequency 178 MHz. Figure|6l shows 
Rdr versus the 178 MHz radio lobe power All along the range 
of the 178 MHz radio lobe power, the distribution of Rdr for 
the quasars is strikingly similar to that of the galaxies. There 
is a marginal impression of a decline of Rdr with 178 MHz 
radio lobe power This trend, however, is most likely caused 
by the mathematical dependency "a/x versus x". It disappears 
when plotting Rdr versus Ljr (Fig.|6lbottom). As a quantitative 
check, we divide the luminosity range into two bins, one for 
luminosity below 210'^^ Lq and one above, and calculate for 
each bin the mean Rdr values. They lie in the range of 10^ ''='=°-3 
(TableO. Again, a possible marginal decline of the mean Rdr 
with increasing luminosity is not yet statistically significant and 
we discuss possible reasons further below (Sect. l4.23t . 

With regard to the "matching strategy" applied by van 
Bemmel et al. (2000) and Meisenheimer et al. (2001) this 
means that, if we choose a galaxy - quasar pair matching in 
178 MHz radio lobe power, i.e. lying in a vertical bin of Fig.|6j 
they exhibit a similar Rdr, exactly as predicted by the unified 
schemes. Some exceptional cases are discussed below. Also, if 
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we plot Rdr versus redshift, or if we use the MIR power at 15 
/xm (with values given in Table|2} instead of the FIR power at 
70 fim, all along the abscissa the Rdr distributions look quite 
similar for galaxies and quasars. 

As illustrated in Fig.0with values given in Table|3] the PR 2 
NLRGs and quasars from Fig.|3] for which we obtained only 
upper limits for the dust luminosity, populate the same range of 
Rdr as the detected sources from Fig.[2do. Obviously the upper 
limits are so high that their constraints are consistent with our 
conclusions on the unification. Future data with the Spitzer and 
HERSCHEL satellites may consolidate this. 

The range of Rdr values spans a factor of about ten around 
the average Rdr ^ 100 as indicated by the dotted lines in Fig.|6l 
This appears large at first glance, but it becomes actually small 
with regard to the range of more than four decades of power 
In view of the complexity of the objects, with a diversity of 
the strength of the thermal IR dust bump, evolving giant radio 
lobes (Kaiser et al. 1997, Blundell & Rawlings 1999) and vari- 
ous circumgalactic environments (Barthel & Miley 1988), such 
a dispersion is certainly not surprising. Therefore, the strik- 



Fig. 7. Same as Fig.|6l but for the sources from Fig.|3lwith IR 
upper limits only. 



ing similarity of the Rdr distributions for quasars and galax- 
ies provides strong evidence in favour of the basic role of the 
orientation-dependent unification . 

4.2.3. Other influences 

When looking at the radio maps, it becomes clear that the jet 
and lobe phenomena are of course also influenced by the cir- 
cumgalactic environment and by evolution. Therefore it ap- 
pears not appropriate, to consider and apply the three "prin- 
ciples" unification, environment and evolution on an exclusive 
either-or basis, rather the interplay of all three principles will 
help to understand the commonalities and differences of galax- 
ies and quasars. The task now is to find strategies for disentan- 
gling the effects of the three principles. Firstly, we discuss our 
data in view of the receding torus model - this topic may be 
considered as evolution of the torus with luminosity. 

4.2.3.1 Receding torus model 

In this modification of the general orientation-unifying 
model (Lawrence 1991, Hill et al. 1996, review by Simpson 
2003 and references therein) both the galaxies and the quasars 
are basically surrounded by a dusty torus, but the inner wall of 
the torus recedes with increasing luminosity of the central en- 
gine due to dust sublimation (rgub oc LP-^, e.g. Netzer & Laor 
1993). Under the assumption, that the scale height of the dust 
torus is independent of luminosity, the opening angle of the bi- 
cone increases with luminosity. 

When defining those objects with direct view to the broad 
emission Une region simply as quasars and all others as galax- 
ies, then the receding torus model predicts that the quasar frac- 
tion increases with luminosity. In fact, such a transition of the 
quasar fraction from about 20% to higher values is observed at 
a 151 MHz luminosity threshold of about 10^''^ W/Hz (Willott 
et al. 2000). However, in all the luminosity bins above 10^^'^ 
W/Hz the quasar fraction is 40% at the 1-cr level (Willott et al. 
2000). This constancy of the quasar fraction at high luminosi- 
ties is surprising in the frame of a pure receding torus model. It 
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suggests that above a luminosity threshold - possibly around 
10^^'^ W/Hz - the geometry of the receding torus does not 
change significantly and that the bi-cone opening angle has 
reached a maximum value compatible with other constraints. 
All of our FR 2 sources have 178 MHz luminosities higher than 
1026 5 w/Hz (Fig.|6ll, hence lie above this threshold. With re- 
gard to our 3CR data we emphasize two points: 

1) If galaxies and quasars belong to the same parent popu- 
lation, then - irrespective of whether a change of the source 
geometry with luminosity occurs or not - any orientation- 
dependent unification predicts that for a given luminosity bin 
the dust-to-lobe luminosity ratio Rdr should be the same for 
the galaxy and the quasar. If the galaxy had a lower Rdr than 
the quasar, this would point to a weaker central engine in the 
galaxy. In such a case, however, we cannot determine from 
the current data, whether this difference between galaxy and 
quasar is accompanied by different torus properties - and even 
if we could, this information would be only of moderate rel- 
evance with regard to other effects like lobe evolution and 
circumgalactic environment discussed in the next subsection. 
Nevertheless, our 3CR data show a striking similarity of the 
Rdr distributions for the quasars and the galaxies (Fig.|6j and 
any difference or decline with luminosity can be established at 
most at the 1-ct level or less (Table|3}- 

2) If the bi-cone opening angle increases significantly with 
luminosity, then the covering angle of the dust, which absorbs 
the radiation from the central engine, decreases. Then one 
would expect that - relative to another measure of the AGN 
strength - the IR luminosity declines. For example an increase 
of the half-opening angle of the bi-cone from 30° to 60-65° 
results in a decrease of the dusty solid angle seen from the cen- 
tral source by a factor of about two (cos30°/cos65°=2.05). We 
searched for such a relative decline of the IR luminosity in two 
data samples: 

a) If in optically selected quasars the blue continuum from the 
central engine is seen largely free of extinction, then the ra- 
tio of blue to IR luminosity should decline. We examined 
the sample of 64 Palomar Green quasars with good near- 
and mid-IR data and covering four decades in luminosity 
(Fig. 4 and Table 2 in Haas et al. 2003). We find that for all 
bins above Lb=10^^ L0 (corresponding to Mb<-24 mag) 
the mean ratios of Lnir/Lb and Lmir/Lb are clearly con- 
stant at the 1 -cr level with values about 1.1 ± 0.5 . But for the 
low luminosity bins with Lb<10^^ L0 the ratios are higher 
by a factor of about 1.5 and 1.7, respectively, although the 
differences are only marginally significant. 

b) For the radio-loud 3 CR objects one may expect a decline of 
Rdr with luminosity, but our data do not show a clear trend 
at a statistically significant level (Fig.|6l Table|3}- 

Thus we conclude that the effects of the receding torus with 
increasing bi-cone opening angle may be significant at lower 
luminosity levels for Seyfert galaxies and radio galaxies hav- 
ing 178 MHz luminosities below lO^''-^ W/Hz, but that at the 
high luminosity level of quasars the receding torus model in 
its originally proposed rigid form plays only a minor role and 
needs some further refinement. 



One refinement should focus on the torus scale height: For 
the receding torus model the assumption of a constant torus 
scale height is most critical. It may not be justified if the AGN 
phenomenon is accompanied by nuclear starbursts and if the 
star formation rate (SFR) grows with AGN luminosity. Then 
supernovae will blow-up the scale height h of the torus (and 
the disk) by injection of kinetic energy. For example, three- 
dimensional numerical simulations of the interstellar matter lo- 
cated in the combined gravitational potential of both a cen- 
tral black hole and a stellar system yield h oc SFR"-^ oc L°-^ 
(Wada & Norman 2002). If this is the case, then both h and 
rsub scale with L*^^ and the opening angle of the bi-cone does 
not change with luminosity. This refinement would explain (1) 
that the quasar fraction remains constant at 15 1 MHz luminosi- 
ties above lO^'''^ W/Hz, and (2) that the ratios Lnir/Lb and 
Lmir/Lb do not decline at luminosities above Lb=10^^ Lq 
(Mb < -24 mag). 

4.2.3.2 Environment and lobe evolution 

The ISO data suggest, that the FIR and MIR luminosities 
actually yield good estimates for the strength of the AGN. Now, 
under the premise that the unification is basically valid, sensi- 
tive infrared measurements offer to pursue the following strat- 
egy proposed by Barthel & Arnaud (1996): We adopt the IR lu- 
minosity as primary measure for the AGN power and examine 
the strength of the radio lobes. Then deviations in the dust-to- 
lobe luminosity ratio allow, free of orientation bias, to explore 
the additional influence of environment and evolution onto the 
extended radio lobes. In the following, we consider the extreme 
cases of high and low Rdr values in Fig.|6l 

High values of Rdr are expected for sources, which are ei- 
ther rather dust-rich or do not show strong radio lobe activity. 
Examples are the "cool" FR 2 galaxies 3C 321 and 3C 459 both 
with a high FIR contribution possibly enhanced by starbursts. 
Very extreme examples for high Rdr sources are radio-quiet 
quasars lying even above the upper plot range of Fig.|6l 

Low values of Rdr are expected for sources, which are ei- 
ther rather dust-poor or do show exceptional radio lobe activity 
- possibly caused by a special circumgalactic environment or 
by a high evolutionary state of the lobes. For example, 3C405 
(Cygnus A) has exceptionally bright radio lobes (Barthel & 
Arnaud, 1996) and shows a low Rdr < 10, confirming the sug- 
gestion by Barthel & Arnaud, that radio luminosity is not gen- 
erally a good measure of AGN power The ISO data show that 
the dust emission may serve as a better isotropic indicator for 
the AGN strength. Also the galaxy 3C 295 lies at the low end of 
the FIR-to-radio distribution, indicating high efficiency in cre- 
ating the radio lobe luminosity. Note that 3C 295 is located in a 
dense X-ray bright cluster, suggesting the influence of the envi- 
ronment in enhancing the lobe power. Also 3C 380 is peculiar 
in that it shows both a relatively flat radio spectrum probably 
seen pole-on and bright FR 2 lobes. A possible way out of this 
dilemma is that there is an intrinsic bend in the jet (Wilkinson 
et al. 1991). In this case a strong interaction with the circum- 
galactic medium may provide an extraordinarily high radio lu- 
minosity, too. 
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The extreme cases with high and low Rdr show that the 
dispersion of the Rdr distribution in general may be well un- 
derstood in terms of evolution of the radio lobes and by effects 
of the circumgalactic environment. 

If evolutionary and environmental effects are superposed 
on the pure orientation-dependent unification, then it is a priori 
not clear whether the increased dispersion affects galaxies and 
quasars in the same manner, or whether a bias is introduced. A 
detailed investigation of the possible selection effects (Chapter 
3.5.3 in Haas 2001) shows that there is a bias, which limits 
the strategies to check the unification: a radio source, which 
has already been recognised as a quasar due to some criteria, 
will on average always exhibit a higher IR luminosity than a 
source with matching radio luminosity, which has not yet been 
identified as a quasar. Larger data samples are required in order 
to overcome this "strategical bias". 

5. Conclusions 

We have obtained sensitive Infrared SEDs for a large sam- 
ple of 3CR radio galaxies and quasars. In order to check the 
orientation-dependent unification of powerful FR 2 sources we 
considered 17 galaxies and 18 quasars with good IR detec- 
tions. The main conclusions for this sample are: 

1) Both the quasars and the galaxies exhibit a high mid- to 
far-IR luminosity ratio typical for AGNs. 

2) The ratio of IR dust to radio lobe luminosity Rdr is sim- 
ilar for the quasars and the galaxies. Therefore, this re- 
sult provides strong evidence in favour of the orientation- 
dependent unification of luminous (Lir > 10^^ Lq) narrow 
line FR 2 radio galaxies with quasars. 

3) The dispersion of the Rdr distribution might be caused by 
evolution by the radio lobes as well as effects of the cir- 
cumgalactic environment. 

4) The receding torus model in its original formulation is 
not valid at the high luminosity of our sources (Lir <; 
10^^ L0). Therefore we propose as a refinement that the 
scale height of the torus increases also with luminosity due 
to supernovae from starbursts accompanying the AGN. 

Compared with the detections, the upper limits of the remain- 
ing 20 galaxies and 15 quasars are so high that they do not 
provide any deep constraints or conflicts with the conclusions 
above. 
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Table 3. Average values of log (Rdr) for galaxies and quasars 
shown in Figs.|6land0 The means ± standard deviations are 
given for two bins in luminosity of Lir and Pitsmhz- N de- 
notes the number of sources in the bins. The averages are also 
calculated excluding extreme outliers. 
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" excluding the extremes 3C 321 and 3C405 

' excluding the extreme 3C 459 

" excluding the extremes 3C 321 and 3C459 

'^ excluding the extreme 3C 405 

"^ excluding the extreme 3C 371 (BL Lac type, outside Figure ranges) 

■' excluding the extreme 3C 380 
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Table 1. Measured flux densities in mJy as a function of wavelength in /im. Numbers in bold are detections above the 3-cr level, 
the uncertainties are 10-30%. Other numbers represent 3-(t upper limits. The apertures used for each filter are listed underneath 
the wavelength; for 1 20-200 /im the apertures were 90" in the case of maps and 180" in the case of chopped photometry, "m" 
and "c" denote mapped/chopped observations. The previous publications of ISOPHOT flux values refer to: A = Andreani et 
al. (2002), vB = van Bemmel et al. (2000), F = Fanti et al. (2000), HI = Haas et al. (1998), H2 = Haas et al. (2003), M = 
Meisenheimer et al. (2001), P = Polletta et al. (2000). In addition to our photometry, we list also the SCUBA F45o^m and Fgso^m 
literature values shown in the SED plots: They are taken from A = Archibald et al. (2001), B = Best et al. (1998), H = Hughes 
et al. (1993, at 800^01), P = Polletta et al. (2000), R = Robson et al. (1998), vB = van Bemmel et al. (2000), W = Willott et al. 
(2002). 
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Table 1. continued. 
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Table 2. Types, redshifts and luminosity distances, restframe luminosities, dust temperatures and masses. The luminosity distance 
Dl is calculated adopting a A cosmology with Ho = 71 kms^^ Mpc^^, r^mattor = 0.27 and f^A = 0.73. The luminosities Lopt 
(0.5-1 /im), Lnir (3-10 /im), Lmir (10-40 fim) and LpiR (40-1000 /im) are computed in the restframes of the objects from the 
"envelopes" of the SEDs. The luminosities give the total emission from host galaxy and the core. The IR luminosities largely 
represent the thermal dust luminosities, except for some sources (like 3C371 or 3C380), which do not exhibit any prominent 
thermal bump and are excluded from the discussion. The thermal dust emission power i^P,y has been computed at restframe 
wavelengths 15 and 70 /im, respectively. The 178 MHz radio lobe power is computed from the observed 178 MHz fluxes, 
applying a K-correction with a spectral index a = 0.5. The choice of a has only little impact on the resulting lobe power and the 
relations discussed. Rdr is the ratio of dust power to radio lobe power at 15 and 70 /im, respectively. The values of Tp refer to the 
peak blackbody shown in the SED plots, the uncertainty is 5-lOK and can raise to 50% in cases of poor FIR upper limits. The 
dust mass Mdust is estimated according to Formula 6 in Stickelet al. (2000) following Hildebrand( 1983), with uncertainties up to 
an order of magnitude. A "^" marks highly uncertain Lmir/Lfir, Tdust and Mdust values due to poor dust bump determinations 
or upper limit constraints, and these sources are not used in the discussion. 
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12.49 


13.09 


12.60 


3.03 


39.38 


38.94 


37.32 


115 


42 


85 


7.08 


249.1 


QSR L 
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35.91 


453 


368 


30 


8.53 


274.0 


NLRG FR 1 
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26 
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8 


<26 


40 


7.94 


298.0 


QSR FR2 


1.436 


10427 
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10.56 
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0.52 
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36 
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40 


7.07 
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2.51 
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55 


29 


70 


7.22 
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QSR FR2 CSS 


1.574 


11690 


11.62 


< 12.86 


13.56 


13.38 


1.51 


39.80 
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37.76 
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134 


55 


8.64 


321.0 


NLRG FR2 


0.096 


435 


10.60 


10.60 


11.33 


11.36 


0.92 


37.45 


37.86 


34.76 


497 


1258 


35 


7.49 


323.1 


QSR FR2 


0.264 


1329 


11.42 


11.62 


11.53 


10.88 


4.49 


38.03 


37.38 


35.61 


261 


59 


60 


5.98 


324.0 


NLRG FR2 


1.206 


8385 


11.25 


12.12 


< 13.00 


<13.01 


~0.98 


<39.16 


<39.56 


37.55 


<41 


< 103 


~50 


~8.42 



* The types are abbreviated as: 

QSR = quasar, BLRG = broad line radio galaxy, NLRG = narrow line radio galaxy, LERG = low-excitation radio galaxy, 

C = core dominated, L = lobe dominated, FR 1, FR2 = Fanaroff-Riley classes 1 and 2, 

CSS = compact steep spectrum source, OVV = optically violent variable. 
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Table 2. continued 



3CR 




type' 


redshift 


Dl 


Lopt 


Lnir 


Lmir 


Lfir 


Lmir/ 


P(15 /im) P(70 ^m) P(178 MHz) 


Rdr 


Rdr 


Td 


Md 










Mpc 


log[L0] 


log[L0] 


log[Lo] 


log[L0] 


Lfir 


log[W] 


log[W] 


log[W] 


15 ^m 


70 ^m 


[K] 


logM0 


325.0 


QSR 


L 


0.860 


5498 


11.43 


<13.10 


13.07 


12.42 


4.54 


39.54 


38.84 


37.14 


252 


50 


75 


7.14 


330.0 


NLRG 


FR2 


0.550 


3168 


11.08 


11.29 


<11.54 


< 11.62 


~0.83 


<37.86 


<38.20 


36.87 


< 10 


<21 


~40 


~7.57 


332.0 


BLRG 


FR2 


0.152 


716 


10.87 


10.70 


10.78 


< 10.80 


~0.93 


37.14 


<37.38 


35.05 


122 


<211 


~40 


~6.75 


334.0 


QSR 


L 


0.555 


3203 


11.54 


< 12.03 


12.20 


11.68 


3.30 


38.58 


38.11 


36.47 


127 


44 


75 


6.39 


336.0 


QSR 


L 


0.927 


6036 


11.79 


11.93 


< 12.64 


< 12.75 


~0.77 


<38.83 


<39.32 


37.09 


<55 


< 168 


~40 


~8.69 


337.0 


NLRG 


PR 2 


0.635 


3777 


10.66 


<11.60 


< 12.08 


< 12.04 


~ 1.09 


<38.36 


<38.60 


36.66 


<50 


<86 


~40 


~7.96 


340.0 


NLRG 


FR2 


0.775 


4830 


10.77 


<11.75 


< 12.22 


<12.37 


~ 0.70 


<38.51 


<38.94 


36.83 


<48 


< 131 


~40 


~8.32 


343.1 


NLRG 


FR2 CSS 


0.750 


4638 


11.13 


<11.83 


< 12.46 


< 12.89 


~0.37 


<38.65 


<39.48 


36.84 


<64 


<431 


~40 


~8.85 


351.0 


QSR 


L 


0.371 


1975 


11.80 


12.22 


12.27 


11.87 


2.50 


38.68 


38.39 


36.13 


360 


185 


50 


7.30 


352.0 


NLRG 


PR 2 


0.806 


5071 


10.57 


<11.56 


<12.10 


< 12.27 


~0.68 


<38.37 


<38.84 


36.92 


<28 


<84 


~40 


~8.22 


356.0 


NLRG 


PR 2 


1.079 


7296 


11.23 


11.52 


12.02 


12.25 


0.59 


38.28 


38.83 


37.27 


10 


37 


40 


8.20 


368.0 


NLRG 


PR 2 


1.131 


7738 


11.35 


< 12.39 


12.92 


< 12.82 


~ 1.26 


39.21 


<39.36 


37.41 


63 


< 89 


~40 


~8.73 


371.0 


BLLac 




0.051 


224 


10.77 


< 10.70 


< 10.63 


< 10.23 


~2.49 


<37.10 


<36.73 


33.59 


<3227 


<1373 


~40 


~6.11 


380.0 


QSR 


C CSS 


0.692 


4199 


11.88 


< 12.28 


<12.33 


< 11.85 


~3.01 


<38.76 


<38.31 


37.46 


<20 


<7 


~65 


~6.81 


381.0' 


NLRG 


PR 2 


0.160 


758 


10.43 


10.91 


11.10 


10.65 


2.82 


37.51 


37.15 


35.34 


147 


64 


40 


6.70 


390.3 


BLRG 


PR 2 


0.056 


246 


10.30 


10.64 


10.65 


10.23 


2.61 


37.13 


36.76 


34.80 


211 


90 


40 


6.13 


401.0 


LERG 


PR 2 


0.201 


976 


10.72 


<10.15 


<10.51 


< 11.00 


~0.33 


<36.52 


<37.56 


35.67 


<7 


<78 


~40 


~6.94 


405.0 


NLRG 


PR 2 


0.056 


246 


9.95 


10.55 


11.38 


11.30 


1.19 


37.64 


37.84 


37.07 


4 


6 


50 


6.78 


427.1 


LERG 


PR 2 


0.572 


3324 


10.59 


< 10.90 


<11.99 


< 12.29 


~0.50 


<37.92 


<38.85 


36.89 


< 11 


<89 


~40 


~8.21 


437.0 


NLRG 


PR 2 


1.480 


10827 


11.15 


< 12.35 


<13.33 


<13.I2 


~ 1.63 


<39.50 


<39.60 


37.76 


<55 


<70 


~60 


~8.21 


441.0 


NLRG 


PR 2 


0.707 


4311 


10.72 


<11.72 


< 12.26 


< 12.42 


~0.68 


<38.52 


<38.99 


36.81 


<51 


< 150 


~40 


~8.36 


445.0 


BLRG 


PR 2 


0.056 


247 


10.26 


10.87 


10.98 


10.50 


3.02 


37.46 


37.03 


34.52 


855 


323 


50 


5.96 


454.3 


FSQ 


C OVV 


0.860 


5498 


12.50 


< 12.92 


13.53 


12.79 


5.59 


39.81 


39.13 


37.06 


566 


118 


100 


7.10 


459.0 


NLRG 


PR 2 


0.220 


1081 


10.89 


11.07 


11.83 


12.24 


0.39 


37.93 


38.81 


35.85 


120 


908 


40 


8.18 


460.0 


NLRG 


PR 2 


0.268 


1353 


10.39 


< 10.15 


< 10.43 


< 10.98 


~0.28 


<36.55 


<37.58 


35.56 


< 10 


< 104 


~40 


~6.94 



^ 3C381 could be classified as BLRG based on a noisy H^ line, which appears to be broader than the [Olll] A5007 line, but with considerable uncertainty (Saunders et al. 1989). A new 
spectrum obtained at the KECK telescope shows that 3C381 should be classified as NLRG (van Bemmel, private communication). 
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Fig. 3. Spectral energy distributions of the remaining sources not used in the discussion, which are either flat spectrum sources 
or have only poorly constrained thermal bumps. The meaning of the lines and labels is as in Fig.Q For the poorly constrained 
cases, however, the blackbody curves and the envelopes are shown here as dotted lines in order to emphasize that they represent 
only a maximum possible dust contribution. For the blackbody, if not sufficiently constrained by the data, we adopted the average 
of Tp w 40 K derived from the good cases shown in Fig.^ The long dashed lines for 3C 138, 3C 273 and 3C 380 indicate the 
extrapolation of the synchrotron emission towards shorter wavelengths. 



